Improved characterization of elastic scattering near a Feshbach resonance in 85 Rb 
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^s^j ■ We report extensions and corrections to the measurement of the Feshbach resonance in 85 Rb 

, ^ cold atom collisions reported earlier [J. L. Roberts et al., Phys. Rev. Lett. 81, 5109 (1998)]. In 

Q addition to a better determination of the position of the resonance peak (154.9(4) gauss) and its 

• width (11.0(4) gauss), improvements in our techniques now allow the measurement of the absolute 

size of the elastic scattering rate. This provides a new measure of the s-wave scattering length 
as a function of magnetic field near the Feshbach resonance and constrains the Rb-Rb interaction 
(N ■ potential. 
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PACS Number(s): 34.50.-s, 03.75.Fi, 05.30.Jp, 32.80.Pj 



The presence of a Feshbach resonance in cold atomic collisions allows the elastic and inelastic collision properties 
to be dramatically altered as an external field is changed In particular, the s-wave scattering length (a), which 
completely characterizes elastic scattering in ultracold collisions, can be tuned over a large range. In 1998, several 
groups reported the observation of external-magnetic- field Feshbach resonances in ultracold collisions in 85 Rb 
and in 23 Na jjj. These tunable interactions have been used to study the stability || and dynamics |(| of 85 Rb BECs 
with attractive interactions. Interpreting the results of these recent studies and other future studies in 85 Rb requires 
• i-H . a more accurate knowledge of a as a function of magnetic field (B) near the Feshbach resonance than that obtained 
in Ref. In order to improve the accuracy of the determination of a versus B, we extended and improved the 
measurements performed in Ref. [||. 

The cold collision Rb-Rb interatomic potential and hence the s-wave scattering length vs. magnetic field is char- 
'— acterized primarily by three parameters: the triplet scattering length, ax', the singlet scattering length, as] and the 
Cq van der Waals coefficient p|||]. In our previous work, we measured two characteristics of the prominent Feshbach 



resonance in Rb: the magnetic field at which the magnitude of the s-wave scattering length was maximum (B peQ fe) 



£f*l . and the magnetic field at which the s-wave scattering length was zero (B zero ). This was done by measuring the 
I/"") ' relative elastic scattering cross section as the magnetic field was changed. This tightly constrained two of the three 
parameters, ax and as, for a given Cq. However, Cq was much less constrained, even when combined with the results 
from another experiment ||]. At the time, the temperature of the trapped atom clouds was too hot to allow the 
determination of the absolute value of the scattering length M . 

We have now been able to improve the measurements made in Ref. Q in two ways. First, the measurement of the 
stability of BECs with attractive atom-atom interactions (negative a) has allowed a much more precise measurement 
of B zero . This more precise determination showed the presence of some small errors in the previous measurement. 
. ' Second, we now have a far greater range of density and temperature of the magnetically trapped atom clouds available, 
and much better control over those quantities. The reproducibility and calibration were also improved. The absolute 
' value of the scattering length could now be accurately determined. This additional constraint allows us to completely 
O characterize a vs. B near the Feshbach resonance. 

The more precise measurement of the position of B zero was a natural outgrowth of our work determining the 
stability of magnetically trapped BECs with attractive atom-atom interactions ||. When large enough attractive 
interactions (defined by Na) are present, the BECs will become unstable and collapse. The collapse criterion is given 
by the relation A^ cr |a| = (3, where N cr is the number of atoms in the condensate and (3 = lAlfim for our magnetic 
trap. In Ref. f§, N cr was determined as a function of magnetic field down to very small a. A simple extrapolation to 
the point where N cr — oo determined the position of B zero . The value of B zero determined in this way is 165.85(5) 
gauss (G), a factor of six more precise than before. This increase in precision occurs both because the stability 
measurement is more sensitive and it avoids temperature-dependent shifts. 

The new determination of B zero disagrees with the previous measurement Q _B^ ero =166.8(3)G by more than the 
stated errors. A small error in the magnetic field calibration in Ref. Q is reponsible for 0.3G of the difference. 
The remainder of the difference is likely due to temperature shifts of the elastic cross section and the temperature- 
dependent asymmetry in the cross section versus magnetic field near B zero JTo[ that were not adequately included in 
Ref. [||. The new magnetic field calibration also shifts the measured position of B pea k to 154.9(4)G. 
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The absolute value of the scattering length was determined at several values of the magnetic field, using "cross- 
dimensional mixing" measurements as in Ref. Q|. In this technique, a cloud of cold atoms, created by evaporative 
cooling as in Ref. jjT| , is perturbed from its equilibrium energy distribution. The time for the cloud to equilibrate by 
clastic collisions is then measured. In the zero-temperature limit, the elastic collision cross section (a) is a — 87ra 2 . 

Once the atom cloud had been cooled sufficiently, the bias magnetic field in the trap was ramped to near B zero to 
drastically reduce the elastic collision rate. By cutting faster with the rf- "knife" used for evaporative cooling than the 
cloud could equilibrate, it was possible to remove energy in just the radial direction of the axially symmetric magnetic 
trap The width of the cloud in the radial direction was decreased until the radial energy had been reduced to 
0.6 of the axial direction, and then the bias magnetic field was ramped adiabatically to a selected value and the 
cloud was allowed to equilibrate for a fixed time. The size, number, and shape of the cloud were measured using 
on-resonance destructive absorption imaging. This procedure was repeated for different equilibration times and the 
rate of the relaxation of the aspect ratio of the clouds (the ratio of the axial width to the radial width) to equilibrium 
was recorded. 

The rate of flow of energy out of any particular direction, axial or radial, is proportional to the difference between 
the average energy in that direction and the equilibrium value multiplied by the factor T = ^(n)a(v) (i.e. ^ i = 

f n(x) 2 d 3 x 

—T{Ei — E avg ), i=x,y,z). In the factor T, (n) is the average density defined by ^ — — , (v) is the root-mean-square 



relative speed y 16 ^ T , f = 87ra 2 is the elastic collision cross section, and k is a numerical factor. The parameter k 
is weakly dependent on the difference between radial and axial energy but is always within 2% of the value 2.50. It is 
calculated using classical transport theory |Q with the assumption that the velocity and density distributions can be 
described in each direction by a single parameter, an effective temperature that is equal to the average energy divided 
by the Boltzmann constant ■ A Monte-Carlo simulation agrees with this classical calculation at the 4% level jl3| . 

The relaxation of the aspect ratio was fit appropriately to determine T and hence a. T is not strictly constant as a 
function of equilibration time since the density and temperature of an atom cloud varied during equilibration due to 
loss, heating, and change in shape. These small (< 5%) variations were included in the fits to the data. The density 
and temperature of the atom cloud was measured directly along with the aspect ratio in the absorption imaging. The 
density was maintained near 2.5 x 10 10 cm -3 and the temperature was kept near 130 nK. For these cold clouds, the 
energy dependence of the average elastic cross section is expected to be less than 5%, and the fits to the data included 
the energy dependence due to the unitarity limit p4| . 

Data were taken at 5=168.0, 169.7, and 251.0 G. The magnetic field was calibrated in the same way as Ref. [||: 
the rf frequency that resonantly drove spin flip transitions for atoms at the center of the cloud was measured and the 
Breit-Rabi equation was used to determine the magnetic field. Colder clouds were used, improving the precision. The 
spread of the magnetic field across the clouds was ~ 0.2G FWHM. Figure 1 shows a set of aspect ratio equilibration 
data and the corresponding fit. Across the range of magnetic fields measured, the equilibration times varied by about 
a factor of 30. The distance from B pea k ensures that the spread in a due to the spread in B across the cloud was not 
significant and that the precise location of B pea k is also not a limiting factor in analyzing the results. 

The uncertainty in the determination of a is dominated by the determination of number and the fit to obtain T. 
The number determination relies on the calibration of our absorption imaging system, which was performed while 
taking into account the frequency and intensity of our probe laser, optical pumping effects, and the small Doppler shift 
arising from scattering multiple photons. We estimate the uncertainty to be 10%, which results in a 5% uncertainty 
in each determination of a. The uncertainty in the fit to obtain T is primarily caused by small oscillations in the 
aspect ratio due to the vertical and horizontal radial directions not being perfectly degenerate. Although this limited 
the accuracy of the fit, the size of this imbalance was varied and no significant shifts in the results were observed. 
Also, the period of the oscillations does not depend on the value of a and so the agreement between data sets of very 
different equilibration times also indicates that this effect is not distorting the results significantly. The measured 
value (in atomic units) of a is -342(10), -97(6), and -60(4) at B = 251.0, 169.7, and 168.0G, respectively The listed 
uncertainty does not include the common uncertainty due to the atom number determination. 

Near the Feshbach resonance, the change in a with B can be well approximated by the relation a = at, g (l — 
fi^flperi ^ [p"5[ . Using this, we can determine a,b g for each measured value of a (Fig. 2). The weighted average of 

all of these data gives ab g = —380(21) where the uncertainty is dominated by the determination of number. 

Our measurement now constrains all three parameters needed to describe the Rb-Rb potential for ultracold colli- 
sions. Using the same methods presented in Ref. 0], the Rb-Rb Born-Oppenheimer potentials are adjusted until the 
theoretically predicted B zero , B pea k, and a^g match the observed values. The calculation of a vs. B is performed using 
a quantum defect method |16|. Matching the Rb-Rb potential to the observed scattering lengths implies ax( 85 Rb) 
= -332 ± 18, a s ( 85 Rb)=3650±g^° @, and C 6 = 4660 ± 20, where all of the values are listed in atomic units @. 
This value of the Cq coefficient is consistent both with an ab initio calculation |19| and a recent photoassociation 
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measurement |20| . 

As a result of improvements in our experimental techniques and the measurement of 85 Rb BEC stability in the 
presence of attractive interactions, we have been able to improve our characterization of the prominent Feshbach 
resonance in 85 Rb. 
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FIG. 1. Relaxation of the aspect ratio to equilibrium. The data shown here were taken at £>=169.7G. Each data point (•) 
represents a single destructive absorption measurement. The solid line shows the fit to the data used to determine the elastic 
scattering cross section and hence the scattering length. The deviation of the data points from the fit around 3-4 seconds is 
due to the oscillations between the two radial directions as described in the main text. 



FIG. 2. The value of af, g implied by several cross-dimensional mixing measurements. This value was determined by dividing 
the s-wave scattering length determined at a particular magnetic field by the factor (1 — Bz g '° B B '"^ ak ). The s-wave scattering 
lengths measured at 251. 0G were divided by an additional factor of 0.99 to account for the deviation of the approximate 
description of a vs. B from the full theory calculation. The error bars shown include the statistical errors plus the fit 
uncertainty due the initial radial energy imbalance. The fit uncertainty is non-statistical and so the central values are clustered 
together more than the error bars would suggest. 
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